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^rl» first pert of this thesis deals with the problem 
of determination of stable region© of operation of a reactor 

l |fo, |jjj^ | 

and second with the represent et ion of certain components of a 
dual cycle boiling reactor poser plant for dynamic studies 

ffh® maximum permissible moderator and feel temperatures, : : 
their time constants, etc. in a nuclear reactor will depend 
upon the choice of the fuel element* lattice structure, etc* )( 
Further , during the operation of a graphite poser reactor the y ■ 
moderator temperature coefficient may vary from -2 ac 10**^ to 

I? \l^r\ | 

15 at 10~ i k /®0* fbtt»0 et the preliminary design stage (in 
view of the uncertainties: in seme p ammeters) the mathematical 
model of a reaotor can be reduced to e nultip&raaeter linear 

| 

syotem. The conventional methods of determination of the 
stability boundaries of a aultiparaaeter linear system ®r® ! 

either too tedious or time c«suat«g or both. It he® been I 

I 

illustrated in this thesis that Idle D-Partition method when 
used in conjunction with a digital computer yield® the stability 
boundaries mush more easily. V”* 

yia the second part of the thesis mathematical models 
for the boiling water reaotor cor® and the secondary steam 
generator for digital simulation are discussed, the mathematical 
model for the reaotor ©or® incorporates the delayed neutron 
effects and the reactivity feedback term® due to the variations 
in the fuel temperature and the voM fraction* 

y The mathematical model for the secondary steam 
g«tl«r*tdr le developed with the aesusptiim tot It sen he , ■ 
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m mm&JMkM-w mL i 

F • Ratio of deviation la -neutron flax to 
steady-state neutron flux. 

X«o “ Steedy-stete xenon- 1 55 concentration fa 
infinite flux. 

X, 0 m Steady-state xenon-l 35 o©»0«®tr»tiQtt to 

steady-state flux divided by Ic* m . fkl ?* — _- 

+ < 5 * fit 

t ,m Deviation in xenon-135 concentration 
divided by , 

I * Deviation to iodine-135 concentration 
divided by . 

©I * Differ ease between statistically weighted 

aeon fuel and coolant inlet temperature® at 
full power, steady-state, ®C. 

Sg * Difference between statistically weighted 

mem moderator and coolant inlet temperature® 
at full power, steady-state, **€U 

* u m Deviation of mean fuel temperature, °0. 

f_ * Deviation of mm moderator temperature, *0, 

Hi 

f| » Fuel time constant, secs. 

fg • Moderator time constant , secs. 

^ of a • m* «ry 1 1 j rtiiln ira ^ mi Ulan 

#s q « £ cgc jjl ior turn zxux, neutr cm®/ cm —sec. 
oC u » Reactivity temperature coefficient of fuel, 

j^®CU 

o<. «* Controller gain per seo. per *0* 

oC m Reactivity coefficient «f »nc a >- 135* 

Hi 

^“x * ' cr oinwe ctica of xenon-135 for 

the rmal n eutron capture, e®*% 



^ l * Radioactive decoy constant of iodine -135, 

' *m f 

mo* . 

> x * Badioeotive decay coast ant of xeaon-t33* 

*» y 
a » 

k « Multifile atlon factor. 

k e » Qhmm in k due to the operation of a control 
rod, 

a * Laplace *8 operator. 



1 ** Are®, aq. ft. i 

°i - C«»«tr.tio» «f i«> ercp of a.!.,*, 
nude 11/ • 

$ * teiperature , °F, 

U « Overall heat transfer coefficient , B.Th.fJ*/W#~®F, 
# 

1 « Heutron generation tine, sees# 

a « Neutron density, n©utron«$/asr\ 
o<r T * Void fraction fioient of reactivity. 

©C* « Fuel temperature coefficient of reactivity, & j 
fi » total fraction of delayed Matrons. 

^ » Fraotioa of delated neutron* of 4** group* 

- Xl S deesy constant of i tis group of delayed 

| 

iflL qjfe» * i 

mm wFwxisi f ii • 

MBtttmxna 

U a* 00l&mi1» # 

ww * ***- flPav ^BpwffPpSr m ♦. '£•?-.«?■ • ... '•.: l 'Vd..' 




SAT m Saturated vapour. 


Co) 


T • Temperate m * «F. 

¥ m Mean temperature, °F. 

W » Haight flow, lbs ,/m®* 

x m Blstame from eat 1*2 nee of ©teoMary beat 
excbmger to boiling point* ft. 

M m las© per foot length, lb®. . 

? m Stew pressure. 

h « Heat transfer coefficient , B.Th.U./eec.- f *f. 
p » Perimeter (in association wit h h), ft. 

^ - tmUm* taat, B.?h.U./lb. 

m>330BI«8 

p «* Primary side* 
s « 3@eoaia.acy aide* 
w * Water. 

8 ** Steem. 
i » Inlet, 
o * Outlet. 


for denoting temperatures three indices notation It 
used. The first shoes the primary or second ary side* the 
second denotes the phase* sod the third denotes whether It 
la input or output <p» entity* Thus* 

Inlet temperature of the primary water* 
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Xa ft nuclear power reactor two typee of instabilities 
am occur, namely (1) fundamental instability, when tin total 
reset or power rites or fall® without any appreciable change in 
the ehape of the power distribution within the acre and (11) 
spatial instability, due to the power rising in mo zone aid 


falling in another oiron though the total roeotor power is constant* 
the fundamental instability am bo eon trolled by means of a 
control rod anywhere in the reaotor whereas the ©outre! of spatial 
instability requires control rods at different points [f] « 
m a graphite power reaotor these instabilities are 
caused by changes in reactivity duo to (i) fuel temperature, 

(ii) mod ©rotor tesperature ami (ill) xenon variation®, fh# fuel 


tonporaturo ©©efficient of reactivity with natural uranium ia 
•2 X to’* 5 ^fc/°C, mainly due to doppler broadening of the 
uranium-238 resonance absorption* the moderator temperature 
coefficient of reactivity may increase from *4 x up bo 1 

♦IS * 10** 5 ^h/*a at high irradiations of 3000 MW days per tame, 
duo to plxttoniiMr*259 produced in the reactor by neutron absorption 
in mm*m*S*f ** mftditM there »r# time lego between change 
In power end oliai«@s in fuel and moderator temperatures [l] . ^ 

Renee ter present stability analysis a point roaster medal 
tawr^rstteg all tease affeete tea baan <*oaan. 

: : llli 




*t m following methods are dtscri bed in the liter store 
fer tlie conventional stability MM lysis of linear systems La t© f] 
1) loutb-Hurwit® criterion, 2} Xy quiet criterion, 3) root locus 
techniques end 4) analog computer solution* then there are 
several variables* the amount of labour involved in applying 
any of the above techniques i# prohibitively large* 

for instance, with m »®log computer the stability 
boundaries «u tw per taster pleas have to be determined as 
follows i all the variables, except two say oC,_# JC*. on ©hose 
Plane the stability boundaries are being plotted are given 
numerical values* The potentiometer corresponding to aC a la 
set to a particular value and the value of cC k corresponding to 
stability boundaries le determined by a trie! and error process* 
Shis is repeated for different values of <<_# this constitutes 

m 

®m set of results corresponding to one ease* this process is 
repeated for as way oases as desired* the stability boundaries 
for one ease by analog computer tabes almost half a day whereas 
by using a digital computer (XBK 7044 ) and IMFarittioa method 
s«e stability boundaries earn be evaluated in 10 sees* (execution 
tine), - 


the equations which describe the performance ♦If %lw 
RM*r Mid at* a* follows [l] > 

^u*u * v A * “V * *o * 0,07 •* <M| 
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1?he modification of the present model have been 
described la the earlier references [l m& 8] . ill the shove 
aquations are linearised hy considering perturb at ions about the 
steady state levels corresponding to full reactor poser, la the 
above model* e mm neutron lifetime of 10 milll-eeoonla hm hem 
rnmmmAf rather than the one mini-second life time of a prompt 
neutron* to take into sow account of the effects of the delayed 
mmmmmBm wm% fNixsMiiMiiii ji^§ #|t W| mm Wg in in# 
mathematical model of the reactor system oen vary* Hence the 
reactor control system designer may like to knot? the variation 
of tte stability bousd striae o n tho plan© and natural 

frequencies of oscillations in me system with the variations in 
the parameter© $ Q> G 1# 0 g * f| and f g . finis the praamt problem 
reduces to tht stability analysis of a multiparous tar ijoeeip system* 

).5 {*% SSili &. '°) « 
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t$li« presence of oil the roots on the loft half plum ensures 
stability , i.e., the Imaginary ®*i® sots m the bound 0x7 for 
the stable and unstable region®, In tlio XMtartltion method the 
imaginary axle is mapped on © too parameter piano of interest, 
so that the possible stable and unstable regions oaa bo identified 
on thie piano* 

with the help of the mathematical model of the gloom 
linear system, tho ohirsotoristic equation is obtained («oo 
Appendix 5*2)* thorn is substituted in tho ehameterietlo 

equation and real and imaginary parts are separately o ousted to 
sera. Set the resulting equation® be 

oC a S t (w) ♦oC^ M ♦ H, (w) • 0 (1.8) 

^ * <£|j C 10 ) ♦ Rg (w) » © (t»$) 

where <£ and cC nre the two parameters on whoso 
piano the stability bouaiario® have to be plotted, them for 
each value of w f and oC k esm he evaluated from the 
following fo rmula e $ 


A * ®| ^2 **^2 




tlifi curve* move along the bcnndary in the direction of ^ -increasing 
and shad • it m Hi® left edge st those points for which 0 an d 
on the right edge for those points for which a< 0. the shaded 
region 1® labelled m stable region in fig®, 1, 3, 5 t 7 sad f of 
this theeie. the only strict method of verifying whether the 
shaded region is a stable region le to apply Booth* e or any ether 
criterion for any convenient point inside this region, fhe rules 
for shading the iwpartition in some special oases are discussed in 
reference® [9 end tOj • 


i.4 


m .Appendix 5*3* the formulae for the evaluation of 
oC a * <<&• ^ and f (for any assumed value of ^ ) are furnished . a 
digital computer .prog rente was written for the above purpose, 

She input to the digits! computer consist® of 0 &$ ® t » 0 g# f| t f g# 
sterling and final values of u> and interval of to # the output 
from the digital computer is obtained in the following format* 

£0 oC « j- A « 

* * (period In hears 1 


These result* have been plotted in fig®. 1 to If* tigs, 1*. 
5* $» T and 9 show tbs stable and unstable regions of operation - 
for different eases* tigs* 2* 4» §# 0 and 10 furnish tht parted* 
of oscillations ®t the different boundaries* Ihstdcmts l l y the 


stability boundaries indicate whether the system is emsditieanlly 
stable or net* am gain margin of the system can ales be ebtalncd 
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Zheorsttoolly spsoMb* » 60 mm taka values from - o> to ♦ » 
In actual practice only © limited region will b© of interest, 
el las# aoat parameter® vary wlHUa oertato Monads, la tfe# present 
ease i oC a can vary to o graphite moderated power resot or frea 
*2 at 1<T* to t5 x to* 5 6 fc /«0* From til® first digital computer 
win, an initial estimate of w we© made for variation ofoC 
between these Halts* Bur tog subse<p«at computer runs « was 
varied around this region* 
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Figs* t and 2 show toe variation of stability boundaries 

end periods of oscillations respectively. Warn tj is reduced 

■ ; • 

Jfep jto 1® rtm -ngni Mgak. 4#%#% a to ggMk •*$?& %ml jfl fc A ‘lit'' iUMltli iHTfc *1 sHij 'tlliljn Jtim Jft,® SM& aMh gt* ssfe'fcb. fill li titffh -rt tftei 

10 «P s®0S# iff®® 20 in# gmB ®l Hi® ttppMr 

•1ABU4V IS m&w% IBif ptiiocl® #1 

o so Illations ere reduced to 2/9 of their earlier values* fines 
tbs results st toe lower boundary ©re not effected# it may be 
oonoluded tost variation of fuel tl»» constent fj bat as sffeet 
to -lower stability boundary* . 




effect of the moderator tim mmtmt® to 2000 mm* 

is to increase periods of oscillations, coapired to those when 



m 500 mm * 



Sine© § t la the difference between fool and coolant 
iolot temperatures, this effect cm he studied bj § % 

suitably* figs. 5 and £ ©how the results for this ease* With 
increasing 0 t> the stable region of operation la i ncre&tiei due 
to the downward shifting of th® lower bounder? « fhe periods of 
oscillations are reduced to 0*016 hoars with Of * £75 o spared 
to 0*030 hoars with Of * 275. The ©oatrolXer gain at the upper 
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stability boundary is not affected 


I ; i I 
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operation i» decreased with increase in by upward shift of 
the lower stability boundary, to upper stability boundary is 
not effected, fig* 8 chows that the period® of oscillations 
mm the lower boundary deoxeee# with tore as lag <* 2 * fh^r ere 
not altered mm the upper boundary, thus affect© only the 
lower atability boundary «M not the upper one* 

*•6.5 jjjiimi im mm ^ 

the effect of variation of stoltoiMi flu* level m 
fto amenta! aui spatial totabilitiee in a reactor have been 
investigated in reference [?] . Only to ftoasental instability, 
can be investigated using a point reactor model similar to the 
present one. The equations of the mthwatical model in Section 
(1*t) can be obtained from to kinetic equations of to model 
in reference [?] by assuming tot to power excursions In both 
regions are ©dual to in phase* Hence it is reasonable to 
expoet that to atability region a to period* of oscillations 
agree with those obtained in reference [?] . which indeed is to 
case* Hence these results are not plotted here. 


1 * 6.6 umLstuamm m m m m m 


A* would be expected to result® are a combination to ’ 
to variation* in 0 t and' §g. figs* 9 to 10 tow toco result*, 
to stable region of operation 1* increased by d toward toft of 
to lower boundary tola to upper stability boundary is not 
affected* to periods of oscillations mm to tour stability 
boundary taw » hit* those at t he uoner boundary hats* 

decreased* 
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In the design of a plant control eastern, in which 
there are several perms* tere that are subject to mrlisbien* the 
advantages of repartition ©ethoi for computation of stability 
boundaries* compared to other approaches is beet illustrated by 
the ease discussed earlier. At the ®m® tine it east be pointed 
out that much of the simplicity of this method is lost if the 
two parameters on whose plane stability bound aria a are being 
plotted appear as nonlinear terms* further information regarding, 
transient response of the system cannot be obtained easily by 
thin method* though an approach has already been made to evaluate 
the transient response using this method in reference [to] • ter 
stability boundary cel eulat ions other approaches like root locus 
technique* analog computer solutions » etc* are rather tedious* 

But they «re useful for the synthesis of the control system* 
taking into account the desired transient response* Hence it any 
be preferable -to use in the prellaimiry design stage the 
©•Partition method, ha the final design stages* when choice of 
most of the perimeters is nearly complete, other approaches c an be 
used for control system synthesis purposes* 



pmp or natural oanvaotion say b« uasd* file street cycle reaotor 
lias Vm disadvantage of not foil awing the turbine load fteswwl* 
to more steam to ©ailed tor, th© re a© tor pressor© i© d©ereaaa§t 
time increasing the ©team void fraction. m& d©ex©a»it$ tin 
r©aatirlt;jr v although m inoreaee ie actually required* In th* 
duel ©yole ay® tea, part of tha ©»a®F fro® 1h# reactor tome 
•teem directly whtoh is fed to the turbine, but part of tin lot 
water prodiieod In the rwaotor Is picked t® a (secondary) a team 
generator* Additional at earn is formed in this generator and ie 
need to supply the turbine, tbs oooied water being returned t® 
the re not or* to increase to turbine demand thus results ill a 


decrease to tenperature of the water and, ®©na**pk«atiy t en 
increase to reaotivity. Because of it® nuaereue advantages » 
meet ef tbs betting mter reaetere are being built with dual-* 
eyele prtotipla*. the 580Kt pleat under eenstrastlto at farupam 
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fh® following amplifying sssiteptiosa are wads to 
obtain a roaotor oor® node! incorporating the aaaaatlal featursa* 

1* 1» tbi axial direction. of th« era®, ainuaaMal addition 

of boat *•■ boos assutisci* Boiling mi non-balling 
length® ax® calculated fro® i&lat aubooollng and hoot 
@da««i in an average etonnal par lb* of water* ?h# 
oroea-aaotion of the avarega choanal i® dsduoad on tha j 
tools of a malt oall gaoaatrar* tflatnslana of the 
avaraga aheawal ara shown in fig* 12* [ 

2(a) For the purpose of beat tranaftr the eternal in ■ 
mmomd to to ooapriaing of two region® *• Cl) tha 
non-boiling and* (11) tto boiling region, An boat 
tranafar ooaffioieats Iron clouding to tto ooo&aut 

- s ' ? ,i ■ J; .. / , * * , > . .A 

in both tha region® have to* mmmi, to to Independent | 
of flow and their values have toan aelaulatad wing 
ataadywstat® to at tramfar data* 


8(b) fto boat transfer along tha axial direction la 

considered to to nagliglbla compared to that in tha 

mkmA *8 aI AfiL ■ 

HPM4 CMnMpto&ilSl* 






epaoe variables 


m mmmmt I that such a simulation represents 
appro s*i»te3y tli© total ears response. 

*Bmm the moS.nl is iiaitea to fee study of im4mmutul 
node o®a motions only* r Ih© kinetic e«gtlw of the senator 
©re m follow C.12] * 

4 S k -/a 55 

*if ™* » ♦ Zw 'X* Q i <*e 

1 i*f 
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r.quatian® (2*1) and (2*2) w ©re linearised using 
perturbation theory around the steady state* file linearise® 
equatlona in terns of perterbatisna from steady-state are 
^TNI ®#&<W I 
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Iot«i The noallse&r term 
equation ( 2.4 }J . 



is retained In 


~ft - /S a> AO i (2*5) 

(t « 1 tO «) 

where m ^ /3 ^ 

w » valises fbr \ i *8 eaft /^ 4 *® far 0-235 fhel 
are given In Table 2*1. 
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Delayed leutven 
Srcmp Wo. 


Delay 

Oonstaiit 


1 


i 


Fraetlmi of ■ 
Total Weuir«s* 

1 


1 

0.0124 

&M1 

2 


0.00141 

5 

0*1110 

0.00127 

4 

0.5010 

0.00255 

$ 

1.1000 

0.0007* 

6 

5.0000 



eepe 1 



Ae oan be seen tmm Fig* 12 t the total channel length 

;' . < . ' > ' ■ 

%b divided into two regions, occupied by { 4 } th® steam vapour, 
(11) the water* bet P* and be the average beet generation 


mtm psr unit length of the fuel red to region* I and 2 
respectively. 

*h© equetioas for the region ( 1 ) ©am bo written as 

follows t 

In any eons latent system of units, per unit length 
of Hie fuel rod tmi for ® snail duration of tine dt sees.* 


The heat energy 
utilised to 
rel slug the fuel 
temperature 


Heat generated 
to the fuel rod 
due to fissions 
to region 1 


teat trans- 
ferred to 
the sledding* 


this oan be symbolically written m 

P f*f Of -ft *f 1 * P < - < f t 1 - *.l> *ft ®fo 1 <*•«> 


Similarly the following equation oen be written for 
the rate of rise of the deciding temperature to region If 

p* *o 0. -Si *.1 - «*f1 ” *o1> *ft °ft1 - 

<*„ - *«*> *o. 0.7) 

for toe seoond region , tbs correspond tog equations 
are as follows • 

P t A, s t -b <*« - **> Aft »*, (*. 8 ) 

P© *© % "If f ©t • ? %2 • f ©f 5 **© %«a m 

(f al - T„) A ot «.») 

tte wrarag* toaporatwe of the water latte am* 

, f . ... ■ ; * * , ’ • 





' ® e ® r® o o 

Tht Talus e far Tarieue heat transfer coefficient e have 
Imw» ealsuluted from the ateaiy-stete heat traafwr data. 




In «ay typieal ohanatl, there *AU te tow region*, 


am ©coupled a&iisly by mater and the ath«r' by met atM* !*«• 



boiling height (» 0 in fig* 13} correspond* to the region in 
eUlih only eeaetbl* heat in added to the xaeming s-ubeooled 
coolant at the Aaa»l bottom. At a • I f » the mrtmb 
baowti* sotureteu . the remainder #f th* obonatl 1© that in 
tfctfill boiling tehee »Um sal la Galled the belling heiaht 


$&m mvtMtso&lma boiling a®y of mm m ccotir in I but © m b* 
neglected for ell pro© ft teal purposes# fha ratios ©aft 

fL*/& may be evaluate a from the ratio of sensible beet added 

M 

to total beat added in the channel# 

Aeaunung sinusoid *1 heat addition In the axial 
direction of the channel* i msLimm. heat addition ret® per 
foot length of the channel will ooour at the ©enter of the 
channel, lef erring to Fig. 15 

% ** _ ' 

j |, g| ftipiltir Sinnnifu. u. nnmnewirnMi ujiirniin. » SMIQ © 

***' Jf jg j( 


Sins© © 8 1 ® the heat aided per lb# ©f water la the 


wbooolei region 


H 


• j \ m **** 


or 


8 


TV 




Substituting for and simplify la# 


• l/t 1 «* Ooe ("* WqfH) 



2,a.3 Mmjmmjs-Mm m j m jmmm* 


equation 


Void fraction &% »s^ m® tics* is defined by tb i 

* - Maputo „„„.„ 

Total volume ©f ate®© water mixture 

in a flow abates es per tho development la ref erenow [Ji 3 


oC W 




41 A- X ’4j| 

* 'iuppppoiip ji j |f| 

g 


(2*15) 


where S * . Slip * Bat la ©f the velocity of the vapour to 

that of the fluid. 

To flag the average void fraction, the total toiling 
length of ft bn«a i divided into a number of aaall segments 

(20 in the present ease). The stems quality footer la 
calculated at each of these sections m& using the relationship 

a© 'AiA ftAg .,***■ '«fe ■* -[rtir Tthmi JP #% ^ JOtk vMHkr«iii m/i* nUitlilft in r sSko ate Atarjahk. wt- ^ dasidM k-. 
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evaluated. The total area under the void f rest Ion serve le 
found by applying trapezoidal rule. The average veld front ion 
0f®y in# snisr# Jtingia ©t lit© ciistti is xotuw ®jr nw 

area by the height of the channel. - 

louatloao ( 2 * 3)9 (2*5) and (2.10) to (2*15) 

constitute the asttesntieal model ©f the reenter cere. 
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in reactivity, there is a rapid rise in the neutron density. 
5h® peak Talas occurs at about 0*2 seas* after which neutron 


density starts aeoressing* fh® average fuel tenpomtnro «3ao 
Increases* a ret flowly md thm rapidly until it reaches a 
mm&mm of 0*B *F after shout 0*7 sees* Mar due to negative 

I 

reactivity feedback fro* increased fuel temperature and void®, 
the neutron density start® decreasing* for a small tine the 
neutron density goes slightly below the steady state level msfi . 
after 2 secs* assumes an almost steady value* the new steady 
state Is shout 0.03 percent higher than the old steady-state 
value* fhue it can he eeen that even without soy external control 

system the boiling water reactor has got self Halting characteristic 

. 
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Pin. 13 shows the variations in neutron densitv 
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average fuel tonporetoro with time, when there Is a step change 





ill inlet subooollag , th© ease plotted here iB obtained when 
there is a stop change of ♦ 2# of mb® oollng corresponding to 
a step change of approximately «O.04f in eoolant inlet 
temperature, She general nature of the response is similar to - 
the previous ease* fhe steady state neutron density is around 
0.0®*f higher than the fall lead vain®. With an increase in load* 
the primary fluid in the secondary heat exchanger loses mere 
sensible heat* Sonseeuently there is an increase in the ears 
inlet subeooliag with increase in load, fhe present results 
indicate that the steady state neutron density (and hence the 
reactor thermal output) increases with increase in Inlet 
subcooling. ' lienee a dual cycle boiling water roaster plant 
hoc got the load following characteristic . • 



a dual eyde boiling water reactor power plant, 
secondary steam generator is am important element. Hence its 
repreaemtat ion is discussed in this section* 


3.1 mmmm mEimim* 


fhe following assumptions Imre been made in formulating 
the mathematical model of a secondary steam generator. 

(i) The secondary ate ms generator om be represented by 
an equivalent counter flow beat exchanger » 

(ii) ffee lOaetio and potential energies «f primary $m& 
secondary liquid - 9 are amall and can be neglected. 

(lit) 2he conduction of beat in the d ireotion of flow ie 
negligible in both the liquids and the natal. 

(iv) The thermal conductivity of the metal in the direction 
perpendicular to the flow is amah that the thermal 
resistance of the metal wall is negligible compared to 
the film resistances. 

(v) Perfect mixing occurs in the liquids each that Idle 

temperature of the lieu id at any section perpendicular 
to the flow is constant, 

(vl) these is no pressure drop in the evaporating region. 

{▼11} the hast transfer coefficient between water and metal 

M — — l — ■ ■— —■ * m.' — — ~ — *. At tttt ( nn in fil if ABW ..mb 

wO \Wkm8 STA@W / 

(tttl) the relationship between tm®, t flux m& the superheat 

QL * fflAw) den he correlated by the expression 
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H*&% utilised Heat transferred 
in raising the * to the metal by 
metal temperature primary U*uid 


Heat transferred 
by the metal wall 
to the secondary 
liouid 


°» \ h ** * *, < V 0,8 *p <Va - V 


- k i <*«* - W 


(5.5) 


Considering unit length along the heat exchanger (on 
the secondary water side) during a small interval of time *dt* 
secs. 


Heat utilised to 
raise the tempe- 
rature of the 
secondary water 


Heat carried 
way by the 
secondary 
water 


Heat transferred 
to the secondary 
water by the 
am tel 


<L f . 
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Based mi the procedure for averaging the temperatures 
in a reg ion explained In Appendix 5.4 
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% similar reasoning, the editions for the evaporates 
nvn 0© ^riiwa Hfl) 




Ittitftfcns (3*2) to < 3.1 t) nr© tfee perfan»&ee ©auction© 
d@8oribiog the ©©eondary boot ©xotanger* for digital simulation, 
ti»e© equation© w©» lineariaea at given in Append lx 5.5. flit 
linearise© equations are given oeiow 

<>.u) 

(5.1S) 
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« A ^*1 “ a 9 “ *10* T m1 + *11 * p bk * *12* ®jw (3-14) 

A %, - * t3 ^T m1 - 4 h * ? aw - * 15 *. * (3.15): 

, ! 

i 

* f pw2 “ A l6 A *pw0 * A t7* f m2 4 A 18* V 4 A 19 a * i3 * U} 

AT pwo * A 20 a %ws 4 A 21 a5 ®2 4 A 22 a W j>w 4 A 23 AX (5*17) 

H A ^m 2 * A 24^ f pw 2 * A 25 Af m2 4 A 26 * ? sw 

+ * 27 ^* 3 ,^ - A je * W p* “ *29 a W bw (3.18) 

** “ * I J0‘^2 + *31 * P gw + *32 * T mi * *33*®«» (3.19) 


3.5 mmmsm SZ ME cOEmcnares or im tjumumm m ?he 

XJjr.iHMtXgBP ECUMIMK5 . 

ft&m the design Amts of the secondly at®®® genera tor t 
following qu antitie® can he obtained, 

(1) Primary water inlet temperature. 

(U) The primary water flow rate, 

(Hi) Secondary water inlet temperature. 

(iv) Secondary water flow rate,. 

(▼) freeenre in the evaporating region* 

Cvi) Hquivalent length of the eountcrflew heat exchanger. 

' Ceil) **» ratio, of the helling to nenboiling length® at 
, full load* . ... : 

{fill) ' Steen tuality factor for the aecondary aide. 


(lx) Tbs ratio® of the heat transfer coefficients fro* the 
primary water to the metal <U t ) ant metal to secondary 
water (tJ^) in both the economising region and tha 
evaporating region* 

One® the steady state values are obta Ined* the 
coefficients appearing in the linearised equations can be 

obt tinea using me formulae given in .Appendix 5*5* 

3*s i mmm m smm » 


tesponses of the required variables, i.e* f the 
•seeadaxy steam generated and me temperature of the primary 
water at the outlet ete* have been obtained for various step 
disturbances* Il|«, (It) thrcugh (21) show the time variation 
in these fUMMiis, fo aimak for the validity of the results 
obtained | steady, state values of 1 h« variables are cheeked for 
consistency as. far as conservation of heat energy is concerned* 
Isf erring te Table 5*1 following observations can be made* 

It is evident from the Table 5*1 that for first 




Change 

laa 

f 

pw® 

•f 


35© 

1 

Ills #/ | 


Hwl# # i 

l 

-10 

il I..J J.O.0I1U. JJPllJ-.IJ L.IJJ JBJ| 

%%&*/ 

«q 

ilfi< 

10 ®f 

i ' :V .. . : 

l - 

f 

-1© 


[ 

$MM$ # 





















From equation (3,6) 

^•w® * *2 f SAf * f swi 

In simplifying equation (3*4), it i® assumed that 

fi; f ®wi ** "k* 1 * s 3if «“* f m t ere slowly varying, When 
ti»r@ i® a step change in there is s step change i n 

8€e ewttoa (3*6) * So the sasurapiion 4y f^ a i® 

negligible is not valid » This explain® the discrepancy in 
result® to Table 3*1 tor stop disturbance to .However, 

©noli step change® are not very common to practise. If this to 
not the ease then the mathematical model sen he easily modified 
m follow®. 

from equation (3*4) 

«* “* h *m» * V* <**»« ~ W> 

0 

* *1 p o ^ W sw* ** ^owf ^ 

■ Tn the above equation the term has to he 

considered. By differentiation with respset to equation (3*6) 
yields* 

S?*«w2 - "2 It 2 »a* * b 2 ft *«t 



Substituting into equation (3.4)* the following 
equation is obtained . 

fc ’“ * ** ** h * "■**«> 

- k i "» <’ f s ») 0 * 8 <\* - »„*> 

fhn® issteM of oaXenlatiag * fro® equation (3*10), 
at is calculated using the above equation* Since 

Sf f ®m " fp f 8Af fc Pt lt em u «^lust®4* eon be 

evaluated If variation lie known fro® the dynamics of 

the other system eomponcnts. 


4 . aaffiBaagg 


3 m 

3. 

4* 


U Haatry, V«li. and Xf ran, 1*1#, "Analogue Ifcvestigatio® &£ 
the Stability- of Graphite Power l®act©r% Proa. I.E.E., 
Vd.110, Bo. 4, April 1963, p* 812. 

Hoffman, 0 * 11 # , "Boir to check linear ays ten stability* • 
Control Itogineering, August 1964 , p. 75 . 

Ibid, part 11, j?©fe. *65, p»84* 

Ibid, part III, Jam * 65 , p. 71 . 

5* Ibid, part If, October *63, p*89« 

6. Ibid, part Y, Marsh *66, p* 81* 

7* Bateses, 2*8#, "Extending South's criterion to check 
ti» relative stability of linear systems* , Journal of 
the institution of Engineers '(India), Yd. 47, Be. 2 M, 
parte October &ml. December If 66 , p. 16, 

8 . Oaebry, Y. 8 , end Deep, G. 8 ., * Stability boundaries variation 
due to delayed neutron . representation in a graphite power 
react©**#. Journal of the ln@titut.lcm of Knglnsers (India), 
Yol. mt $ So* It, pt* 8 X 1 6 , August 1966, p* 571* 

9* Aiaerman, M*A#, ‘’theory of Automatic Gonlrd 8 , Pergemon 
Bress, 1963# 

10# Bearer, M»f* t M Hatred unties to the dyasnies of Aut omatic 
Regulating of ttectrioal Machines** . Suit® rwor the, leaden. 


11 


Xenshroa, 8*1. and Higgins, felt, "2Kree«pe*iti«i analysis 
«£ Antons tic Control Syst on«i H® print Bo. 4#5 of the 
i^imering Sxperiaeat Static© at fee ttoiversife- of 
Wisconsin (0 *»’:*&•) 








1 x 10 J to a * 10 14 mmtatmq/m* 


9mm equation Cl. 5) to metim l.f 


of I in tquatlo* Cl 




■ MM 


<*. * 


Ci-x 0 ) 

f m -ri l l.. I 

3 Xi 


f^l m 




^ jj 


f ■ * 
I> 


X i <’-*<,> 


Etob o Wat ion ( 1 ,t) substituting for <C U , * >f * , 

jpua v 

*™WW^"w’ 'wPPtgjJgiL 

W 


t X 


lli . V* * 0.027 o 1 (»-t^B)> 8 ,^ g 

n nnumijn m|> tMMwWiii* n» — » «*M* " mi+ m t mmmmu* t un t f& mM mnn m m m TTirffiTirti i ffi t OTiM i ^ 


Cl ♦ f f *) (i*f 8 ®) Cl+® 5 *)(l 4 f 4 ») »(i+f 1 ») 

» o#o ?7 »? 

1**1 A| * *2i tor 4 5 

h m ®z 


Ay • • 0*027 


% * -®| 


4 , - - 0.07 

equation m yields* 



*1 [• t (* 2 »Tj+* 4 )«* ♦ (* g *j«* a *^+* j» + )»* ♦t 8 » J T 4 » 4 J 
*V«[ ** <WV* 8 +(* 1 * J +* 1 I 4 -rt,I 4 )« 5 ♦* J I J * 4 ,*j 
♦*> [•+ (VW®* ♦(l' 1 I 8 -f 1 * 5 -S a I 5 )»5 -I,T a * 5 •♦] 

♦VkC u < VW* + <www* a 3*4' 5 J 


>4 J [* 8 + <*i*VV*« ) * 5 ♦< I lV Sf lV T lV S aW4 4, J*4 > i 


*f|W *l*aV *1*J*4 + W* ),S 


» 0* 

5.5 ITAKUflO# Of/ . <* frl * A® ft 

8ufc»titutiB« • « Jw in #i® wot Ion 

«twi*#<t 4 o •action 5.2 ana real as* iw^txuoey fart® 

separately to*®*© ytald® ttie following too twaUauu 

oc a s, C=o) 4oC k cs f » * i t (w) - 0 

of. C 03 ) ♦ <^ k ( w ) ’&> M • o 

4w * m* m m- 


m m m m 
* I *1*3*4 


] 



•l 0") - -4, (*, ♦ Tj * * 4 ) w* ♦ A, *, *J * 4 J 


wllMptl 


®2 ^ m *2 ^ **■%> ♦ $|T| + )<*>•* 

% C^) • * 4 <fg ♦ f 5 ♦ » 4 >W . ^ f 2 * f^o3 

*1 ( w ) » Ij to 5 •Bj ♦ l| w 

( B| to Bg are as definei below ) 

B| * 4 | ♦ 

*a * *9 ♦ *3 (f^fg-tg) ♦ A|(f s 4f^4f 4 > 

% - 4 5 (f 1 ^ 2+ T 3 4l’ 4 >A 3 (f t f 2 -U t f r T 2 f 5 ) 

♦ A| (fgf^fgf^f^) 

. B 4 • MVa^lV^VVl^fV^V 

■ + A|® 2 * 3® 4 

*3 " *3 (, lVj + *1*2*4 * *1*3*4 * *2*3*4 > 

% * A j *|V |®4 

for a giwa value of <0 § $ %t ® |f n n Sg f ^ f and fig 
am Be evelaiatet tmm the above fo melee* thm a * or^ ocg 
antt f tail be ev%l«N§i# — Imm* 

* • % « 8 - h % :.' ■' : 



?arioi f 


* " " f feg p fcwst 


for a typical 0 interflow teat exofcariga*, when tfe#re ■ 

1.8 no ohpMgQ of phase, following Si ffirtallal aenatlou© iaaorlh# 
the teat trm&tw* 


% V is ? 'pw * ~ f 1 <V "V • 
°w *W fe 2 aw * (f * “W • 

rrm tte above two e«matl©na, 

■ O „ hJm I u 2 V 

♦ % 


0.4.1) 

0.4.2) 


0.4.5) 




■-fc* 


vfcv* °* %» fe T «r ■ 0 


<S.4.4> 


XhthStetisg tte aboro two ©ousationa with reapeot t# 
x ate nstng tte vaml boundary ©esaiUon©* 


c * 
11 


pw pw 


#*U 9. 

w 


2 • 
MB ■ 


' 

?,. , : 




■ ■ • ■ 


■ ’ 't'lif ' i 

• •«. . 


CL w 2 a H3L W. 

w JP© pit 9 


f 


■ i ■ ■>■ 


' w& 

(I.4.SI 


Fro* (5.4*5) 


* P* ^pwi 85 awi? (5*4*6) 
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l C( 25) »PL£*t{25 I *ttTS(25 ) »AR(4) »$Rt 4) »CN(4| 

COMMON N»H>T»TO»rOASH»Y.AR*AR0CF*OR 


ISTRT 


A0&FiX)«A8$tX) 

Aft(l)«*«9 
AR!2|«®*2920f322 
A*< * >»0 # I 70T104ACG0 1 
AR(4|«0*IA*6**6T 

MUMItl 

iJt<2f**»0 

cnui«*«ft 

Cft< 2 >*0*29289922 

€2lB»0*f29**989 

C2*E-0.74»801*5B~0* 

V' . 

CONTINUE 
ISTRT • ' 1 ' 





e%n f% f% f% 


C 

2 


U 

201 


2y2 


207 

206 

204 


12 

602 


5 

2220 

222 


101 

lot 

104 


print*©,* 

PRINT30*T *H*F I WALT* HEX IT 
>*INT30»CY<I)»2«1»N) 

PRINT3o*A*A1,A2,A3 

PRINT30,HLEAST 

PRINT30*(WTSt n*I«l,N) 

TQ»T 

CONTINUE 

CALL DERIVE! Y»D£RINT) 

CALL START IPLESCJ 
C25 1 1*0*3 75*H 
SIGMA«0*0 
SUM*Q«o 

00 202 1*1, N 

WTSI«MTS(I) 

SIOMA»SIGMA<f A6$p C ¥ € 14-75 i )#*fT3I 
ERROR* YOASHfii YOASH! 14-25) 

SUH.SUNWTSt.MWF.ru* 75|-y( 1 1-CMI j.^rror^rror^rror 

SUM- 3 UM.a., 62963 + 

IF I SIGMA— 1*0)206*206*207 

SUM-SUM/SIGMA 

IFISUM— A)5*5*2U4 

T*T— H— H— H 

H*0*5*H 

IF!H-HLEASTI 12*201*201 
PRINT 80 
PRINT 30 *T 

PRINT 30*! Y! X ) *I»1*N) 

GO TO 1 
T * T -H-H-H 
PRINT 320.M 
DO 6 I * 1*51*23 

NEXT STATEHENt* 3 T*oCF£I(d*s***pof*the**co«p*lER* **•«**•* 



T «7*H 

— *** 

* — * 

C21A * 1 *3333333#H 
C21D * M4-H4M 
00 104 !*l»M 

$UM*Y0A»ftf IM-Y0A5H! I* 50) 

pui * vc i i » 

fcM! I i*p( n*C21U*PLESCII 1 

T»T4-M 

CALL DERIVE! EM *EM0A5Mi 
IRROR*0»0 




NEXT T*0 STATEMENTS OCPEmd 

JJJf 

CALL OUT I HEX I T * Y 1 7l» 1 1 

IPiriNALT-T»$9*99*i8e 

IF{ERROR~AHU**m»103 

H*H+M 

DO 113 I*1»N 
Yf D'YCI+TSI 
CALL STARTtPLESCl 

GO TO S 

IF{ ERROR—A3) li3*124»l2A 
T*T-H 

IP <H - HLEAST ) 12* 12§*1»« 


upon the compiler 
****•■***»****#*##• 







CALL DERIVEfYf I+l I tCAY) 

IF((U~U + U- 7Si*U~ SOniA.lStlA 

DO 16 IJ»l»N 

JJ«*I J4-I-50 

YDASH( JJ»*CAYIIJI 

CONTINUE 

IF f J— U/2*2I5*3 »5 

T-T+HALF 

CONTINUE 

DO 8 I 1*1 *N 

QTY*AY*{CAYi I I 1-8E£*0( 11)1 
I II*I+I | 

YU1I )*YI III )^OTY # H 

Of II )»ail I)+3.0»QTY-S£A*CAYCI 1 1 • 

CONTINUE 



*** 

THE 

*** 

CAL 

DO 

PLE 

PET 















4, TG*TO+HEXIT ■ 

PRINT 20»T»YU)»TFAVG»tfOlDtHZEROftY<tt} »Y(10* ,Y(9)tYUll 
X RETURN 

20 FORMAT !8X»F10*6*4£16*4»4£1O*3/} 

END 

tIBFTC SUBS MODECK»NOPRNT 

SUBROUTINE DERIVE! Y*YDA SHI 

€ 

DIMENSION SPACE! 213 1 j 

C 

DIMENSION Y<1O0>» YDASH(75) I 

COMMON H • H > T » SPACE » ISTRT | 

C 1 

COMMON ALNTH » MIGHT*HZ£RQ*FLOtf*PAI *QS*C0NST*TSAT*TSW*TC1.TC2*A2 ) | 

1 A4*F AC1 *FAC2* VOID 

COMMON TFAVG I 

C 

AT ANF t X I * ATAN ( X 1 
C 

1 FORMAT (//0X. 8H TIME »16H NEUTRON DENSITY* 16H AVG FUEL TEMP • 

1 SX*6H V0ID*1UX*6H HZ£RO»8X»5H TF1»3X*5M TF2*5X*5H TC1.5X* 

2 5H TC2I 

2 FORMAT !/2GX*25W REACTIVITY DISTURBANCE **£10*3/1 

III FORMAT ( 20X * 31H PERCENT CHANGE IN SUBCOOLING ■ * £10*3/1 

123 FORMAT! £10*3 ) 

124 FORMAT (20X*6H TwI **ElO*3/2iX*5H QS »*£10*3> 

IF i ISTRT - l I 15 * 14 * 15 

C ■ 

14 FAC1 * 0*66237628 , ' 

FAC 2 * 0*60389151 
C 

QS*42* V 
TWI-510. 

READ 123*CNTRD*OS0 

QTY * OSD * QS / 100. ■% 

QS* OS * QTY 
TWI * TWI * QTY 
DELTA * CNTRD 
HZEftO * 5*7595449 
TCi*620. 

TC2-600. 

AN»2.272?273E+0B 
S5P*6l* 182837 
AX«SSP/AN 
C 

01*21*86867 
02*0*057686028 
83*0*28899514 
84*7*40830330 
85*21*868670 
86*0*038794642 
87*0*29434906 


-ni3% 



f i n nn a ft ft A ft 


88-7.0200S160 

CONST*0.G2U/0.4456 

HIGHT-12.5 

TSAT»544#6 

TS**tTWI+TSAT>/2. 

PAI»4**ATANF(i.) 

FLOW*0*62489132 

ALNTH«649*4 

A2*.06/5>92822 

A4»*O64076664 

V0I01- 0*16635496 

VOID - VOID! 

CALL STAR 

PRINT 2*CNTRD 

PRINT 124 » TWI » QS 

PRINT Uli QSO 

CALL STAR 

PRINT 1 

CALL STAR 

ISTRT * ISTRT + i 


YDASHI i ) *»— 6#45*Y ( 1 J 40*012 4*y ( 2 }4U#0302>*Y ( I l40*lil*Y 14)40*301 
1 *Y i 5 ) +1 . 1*Y <6 ) +3*3*Y ( 7 )410O0*#0ELT A* { AN+Y U > > 

YDASMI2 I ■0»3!1*Y 111 - 0*0124#YI2I 

YDASH(3)«U41*YU> - 0.030S*Y(3) 

YDASHI4 > »1*27*Y< 1 ) - 0*111*Y(4) 

Y0ASHt5»«2*55*YU) * O*301*Y<5> 

YDASHt6»*0.74*YUI - U1*Y<6> 

YDASHI 7 ) »0»27*Y< 1) - 3.3*YI7> 


P0W£R«AK*YU ) 

QMAX* ( POWER/HI GMT )*<PAI/2.) 

Pl*FACl*QMAX 

P2*FAC2*0MAX 


YOASHI 6l*»8l*Pl - B2*(YC8)-Y<9) ) 
YDASHI 9 1 »B3* ( Y C 8 I**Y 1 9 1 1 »84*Y f 9 ) 
YOASMUOI* &5***2 ** BG*IYflOI^vmil 
YOASMUII • B7*«YU0)-YC 




TFAVG * U HI6HT * H2£ft0 * * Y(PJ * H imo * YtlOU / WIGHT 
TC 1 « (2 0. ♦ Y| 9 J 

TC2 « goo. + y(jij 
VO I 00 * VOID- VQID1 

OELTA»f«*f.E«*'04) *TFAVG m O*1*V0I0D 4- CNTKD 

RETURN 

END 

SI8FTC SUB4 N00ECX#NOPRNT 
SUBROUTINE CHANGE 
DIMENSION XI2O),V(20UAREA(2O> 

DIMENSION SPACEI21SJ ' ' 

COMMON H 

COMMON SPACE » ISTRT 

l C ^^S.;*cn^;r R 0' FL0W,P * l ’° S,COKiiT ’ TS * T>r5 “'TC 1 -T«.*. 

COSF I X ) ®CO^l m \ 

£**(** A2 * lIC1 " TSATMIHI6HT-H/ER01 

^ 2 .%^ , r^a TSW,#HZER0 

QT* SUM/ FLOW 

i^o1 RC0S< 1 * ** 2**(QS/QT } } 

H/LRO » MIGHT # A / PA I 

CALCULATION OF NEW VALUES OF FA Cl AND FAC 2 

F4C1 - ( 1# 4 CQSFIAJ > / t PAI * A* 

FAC2 * < 1. * COSFIAJ J / A 


HINCT * ( NIGHT - HZERO J / 20. 

CALCULATION OF THE. HEW VOID FRACTION 

QMAXl « ( QT/HIGHT J * | PAI/2. 1 
START » COSF ( A J 
DO 10 J » l * 20 
AJ ® J 

®, * f AI '* * H ! ER0 * AJ # Hmc * 1 2 hight 

FINAL * COSF f ft $ 

HEAT » I «MAXI*MIGW? /PAI J * * START * FINAL » 

XI JJ * MEAT / ALNTH 



A A 


AX * Xf'JI 

C » M U * AX J /AX ) * CONST * 2, 

10 VfUI « 1. / 11, + C ) 

A«EAU*«mi»HJNCT/2. 

00 30 4-2,20 

30 AREAf4i«m4-l)+V*4H*HINCT/2. 

SUM 1-0 • 

DO AO J*l,20 
40. SUM1-SUM1 ♦ AREA( Jl 

VOID * SUM! t HI6HT 


RETURN I 

END | 

SI8FTC SUBS NODECK, NOFRNT | 

SUBROUTINE STAR 

PR1NT10 I 

PR1NT11 

PRINT10 

10 FORMAT { IX ) 

11 FORMAT ( 1X»116H *»###!«##**»*#*#***###•#*#*«**##***#*#*♦##*#•**#*< 
1 ******##*#*■**# • ****•• *»*••»•*#*»**••#•* *-**#•*«#« •***#*»•**»**«#• 

RETURN 

END 


SENTRY 





0 2 4 6 © JO 12 * 14 16 18X10 

MODERATOR TEMPERATURE COEFRlClEt^m} 6k PER DEG. C 
F!6,2' EFFECT OR VARIATION OF "q ON PERIODS OF OSCILLATIONS. 



CONTROLLER GAINJf- K J «* / Stt/L. 


) 



PERIOD* HOURS 





FIG ,4‘. EFFECT OF VARIATION OF ON PERIODS OP OSCILLATIONS 







PERIOD. HOURS 



FIG.6 EFFECT ; OF VARIATION OF G. ON PERIODS OF OSCILLATIONS 




CONTROLLER ' GAIN k), &k/sEc/c 






PERIOD, HOURS 



OSCILLATIONS 



HIGH PRESSURE 
SATURATE!) STEAM 


STEAM 
S.<Ji AKaTOR 




m\iw m 



T 


j.v— 


TU1BXISE 




LOW PRESSURE 
STEAM 


SECONDARY 
ST? AM 
G Eft x,’ RA TOR 


CQHDKMSER 


RECIRCULATION 

PUMP 

—rtnwtw 

\ 

«S> « *s r^f,% 



LOW PRESSURE 
PEER WATER. 







1,0 AD 


fc W n rwsa^TWBsa 


T"~"~ 


OOH33WHSATB 


high pressure 
PESO WATER 


nt. XX Basic me^to 

: ■• ; "'V; ; - B; : ■ ;: : : ? 

•: 


aBaa 






action and non 















tbs*/5ec 







I - t 

gi;, [tS* 


SIEP CHANGE IN W 
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